Microwave potentials and optimal control for robust quantum gates on an atom chip 
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We propose a two-qubit collisional phase gate that can be implemented with available atom 
chip technology, and present a detailed theoretical analysis of its performance. The gate is based 
on earlier phase gate schemes, but uses a qubit state pair with an experimentally demonstrated, 
very long coherence lifetime. Microwave near-fields play a key role in our implementation as a 
means to realize the state-dependent potentials required for conditional dynamics. Quantum control 
algorithms are used to optimize gate performance. We employ circuit configurations that can be 
built with current fabrication processes, and extensively discuss the impact of technical noise and 
imperfections that characterize an actual atom chip. We find an overall infidelity compatible with 
requirements for fault-tolerant quantum computation. 
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I. INTRODUCTION 



The physical implementation of scalable quantum com- 
puting unquestionably poses huge challenges. Many pro- 
posals exist, most of which in principle satisfy all require- 
ments for scalability 1] ■ One might be tempted to claim 
that their realization is simply a matter of technological 
improvement. Yet, we can fairly say that to date no- 
body knows how to build a scalable quantum computer. 
Understanding and overcoming each of the sources of im- 
perfection, specific to a given physical scenario, becomes 
thus crucial on the way to a possible actual implementa- 
tion. 

Atom chips 0.Q combine many important features of a 
scalable architecture for quantum information processing 
(QIP) The long coherence lifetimes of qubits based 
on hyperfine states of neutral atoms Q , accurate control 
of the coherent evolution of the atoms in tailored microp- 
otentials @, |(J , scalability of the technology through mi- 
crofabrication 0, , which allows the integration of many 
qubits in parallel on the same device while maintaining 
individual addressability, and the exciting perspective of 
interfacing quantum optical qubits with solid-state sys- 
tems for QIP located on the chip surface Q. However, 
the experimental demonstration of a fundamental two- 
qubit quantum gate on an atom chip is an important 
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milestone which still has to be reached. 

A first proposal for a quantum phase gate based on col- 
lisional interactions between atoms on a chip was put for- 
ward in |lOj | . While it demonstrates the working principle 
of such a gate, there are problems which could prevent a 
successful experimental realization: 1) The qubit is en- 
coded in two states with a magnetic-field sensitive energy 
difference, such that it is hard to maintain the qubit co- 
herence over a long time in a noisy experimental environ- 
ment. 2) The fidelity is strongly reduced by wave packet 
distortion due to undesired collisions in some of the qubit 
basis states. 3) An idealized situation with accur ately 
harmonic potentials for the atoms was considered [10j, 
which is hard to realize experimentally, and deviations 
from harmonicity spoil the gate performance ■ 4) Fur- 
thermore, transverse excitations of the atoms were not 
considered in detail. For a successful experimental imple- 
mentation, however, a scheme is needed which allows for 
high fidelity gate operations under realistic conditions. 
A gate infidelity (error rate) below a certain threshold is 
needed in order to allow for a fault-tolerant implementa- 
tion of quantum computing. Depending on error models 
and recovery schemes, estimates of such threshold vary 
from a few 10 -3 for active error correcting cod es [121 up 
to well above 10% for error detection schemes |l3|. 

In this paper, we present a substantially improved ver- 
sion of the phase gate of [lOj for a qubit state pair with an 
experimentally demonstrated coherence lifetime exceed- 
ing 1 s |4| , and give a detailed prescription for its imple- 
mentation. The decoherence rates of the qubits due to 
magnetic-field noise are suppressed by a factor of 10 -3 
compared to |lOf. In our proposal, a key role is played 
by microwave near-field potentials on the atom chip Q, 
which allow to create the required state-selective poten- 
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tials for a successful gate operation with the states con- 
sidered here. Using microwave potentials, the same ro- 
bust qubit states can be used for information processing 
and storage. At the same time, the microwave potentials 
avoid the unwanted collisions limiting the fidelity in the 
original proposal, by appropriately displacing the poten- 
tial minima for qubit states |0) and |1). We simulate the 
gate dynamics in a potential created by a realistic atom 
chip, which we specify in detail, and which can be fabri- 
cated with today's technology. Furthermore, we consider 
many sources of infidelity in detail, such as loss and deco- 
herence effects due to the proximity of the chip surface, 
and we find a total infidelity of the order of a few 1CP 3 . 

II. THE COLLISIONAL PHASE GATE 

A. Qubit states and single qubit rotations 

Quantum information processing with high fidelity is 
only possible if qubit basis states are chosen whose en- 
ergy difference is sufficiently robust against noise in re- 
alistic experimental situations. In particular, technical 
fluctuations of magnetic fields are notorious for limiting 
the coherence lifetime of magnetic-field sensitive qubit 
states of atoms or ions to a few milliseconds 01- O n 
atom chips, magnetic near-field noise due to thermally 
excited currents in the chip wires is an additional funda- 
mental source of decoherence for magnetic field sensitive 
qubit states 0|. To achieve long coherence lifetimes on 
atom chips, it is therefore highly desirable to choose a 
pair of qubit basis states whose energy difference is in- 
sensitive to magnetic field fluctuations. We choose the 
\F = l,m F = -1) = |0) and \F = 2,m F = +1) = |1) 
hyperfine levels of the 55i/2 ground state of 87 Rb atoms 
as qubit basis states. The magnetic moments and the 
corresponding static Zeeman shifts of the two states are 
approximately equal, leading to a strong common mode 
suppression of magnetic field induced decoherence. Fur- 
thermore, both states experience nearly identical trap- 
ping potentials in magnetic traps, thereby avoiding unde- 
sired entanglement between internal and external degrees 
of freedom of the atoms. Coherence lifetimes exceeding 
1 s have been achieved for this state pair with atoms in 
a magnetic trap at a distance of a few microns from the 
chip surface Q. Single-qubit rotations can be easily im- 
plemented by coupling the states |0) and |1) through a 
two-photon microwave-rf transition, as demonstrated in 
1. 

B. Principle of the gate operation 

To complete the toolbox for quantum information pro- 
cessing on atom chips . 16'?] , a fundamental two-qubit 
quantum gate is needed. Each qubit is represented by 
an atom in a superposition of the robust qubit states |0) 
and |1). Our goal is to realize a phase gate, with the 



truth table 

|00) -> |00), 
|01) - |01), 

1 10) -> |io), 

|11) - e**-|ll>, (1) 

for the four two-qubit basis states. 4> g is the gate phase, 
which is to be adjusted to 4> g = n. As proposed in [Tof . 
the phase gate can be implemented by modulating the 
trapping potential state-selectively, such that the two 
atoms interact and pick up a collisional phase shift <p g 
if and only if both are in internal state |1). In the follow- 
ing, we briefly sketch the working principle of the gate, 
highlighting the differences between the present approach 
and [nj. 

In Figure ^ the principle of the gate operation is 
shown. The atoms are placed in a state-dependent po- 
tential 

Cr 1 (r,t) = u c (r) + A(t)-Ui(r), (2) 

which can be split into a common part u G (r) and a state- 
dependent part tij(r), where i — {0, 1} denotes the states 
|0) and |1) and r = (x,y,z). The common part of the 
potential is time-independent, while the state-dependent 
part is modulated with a function A(i), < A(i) < 1, 
during the gate operation. At times t < 0, when the gate 
is in its initial state, we have X(t) = and the atoms are 
subject to u c (r) only, see Fig. [IJl The potential u c (r) 
provides a tight confinement in the transverse dimensions 
y and z, such that the dynamics of the atoms is effectively 
one-dimensional. In the longitudinal dimension x, u c (r) 
is a double well with a sufficiently high barrier to prevent 
tunneling between the wells. Each of the qubit atoms is 
prepared in the motional ground state of one potential 
well. The gate operates during the time < t < T g , 
where r g is the gate time. During this time, A(t) ^ 0, and 
the potential is state-dependent as sketched in Fig. ^p. 
The effect of Uj(r) is twofold: tii(r) removes the barrier 
of the double well for state |1), leaving only a single, 
approximately (but not exactly) harmonic potential well 
in which atoms in state |1) start to oscillate. «o(r) shifts 
the minima of the double well for state |0) further apart in 
the x-direction. The effect of itj(r) on the tight transverse 
confinement is very small. In this way, the truth table 
is implemented: In state |11), both atoms will oscillate 
and collide each time they pass the center of the trap, 
which leads to the desired collisional phase shift of the 
state |11). In states |00), |01), and |10), the atoms do not 
collide, since atoms in state |0) are shifted out of the way 
of the oscillating state |1). When the desired phase shift 
<f> g is accumulated after an integer number of oscillations 
N of the state |1), the gate operation is terminated by 
returning to A(t) = for t > r g , recapturing each atom 
in one of the potential wells of u c (r) (see again Fig.^). 

The gate operation described here is different from the 
original proposal in 10]. There, only the potential for 
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FIG. 1: (Color online) State-selective potential, atomic wave 
functions, and principle of the gate operation, (a) The state- 
independent potential u c (x) along x for t < and t > T g , 
before and after the gate operation, when X(t) = 0. The 
initial state wave function of the two atoms in this potential 
is shown, (b) The state-dependent potential u c (x) + Ui(x) 
(here X(t) = 1) for < t < r 9 , during the gate operation. 
The atomic wave functions after half an oscillation period 
are shown. The state-independent part u c (x) is shown for 
comparison. 

state |1) is switched during the gate operation. This 
leads to unwanted collisions of atoms in state |01), which 
are a major source of infidelity and would require an ad- 
ditional transverse shift in the potential minimum |lfj| . 
Using microwave potentials, we are able to create a state- 
dependent potential in which the state |0) is shifted out 
of the way of the oscillating state |1) in the longitudinal 
dimension. This is essential to avoid the unwanted col- 
lisions and to achieve high fidelity gate performance. A 
second difference is that we choose a smooth modulating 
function X(t) for the state-selective part of the potential 
instead of instantaneous switching as in {lQj . Roughly 
speaking, X(t) is ramped from X(t) = to X(t) = 1 dur- 
ing the first half of each oscillation, while it is ramped 
back to X(t) = with the inverse temporal profile during 
the second half. This smooth ramping allows for much 
better control over the gate dynamics than instantaneous 
switching of the potential. This is necessary to avoid 
strong excitations of the atoms in |0) during the shift of 
the potential well. It also decreases the oscillation am- 
plitude of state |1), thereby further suppressing collisions 
in the states 1 01} and |10). The exact time dependence 
of the function X(t) is determined by an optimal control 
algorithm which optimizes the fidelity of the quantum 
gate, as discussed in Section IvTl 

III. MICROWAVE POTENTIALS 

The state-dependent potential (0) which is needed for 
the gate operation can be realized by a combination of 
static magnetic and microwave fields on the atom chip, 
as briefly described in Q. The concept of microwave 



potentials is similar to the optical potentials created by 
non-resonant laser beams, which have been used with 
enormous success to generate potentials for the manipu- 
lation of ultracold atoms The main difference is that 
spontaneous emission is negligible for microwave transi- 
tions, so that one can use (near) resonant coupling. 

In 87 Rb, microwave potentials derive from magnetic 
dipole transitions with a frequency near ujq/2tt = 
6.835 GHz between the F = 1 and F — 2 hyperfine man- 
ifolds of the ground state (Fig. |2J). The magnetic com- 
ponent of the microwave field couples the \F = l,mi) 
to the \F = 2,7712) sublevels and leads to energy shifts 
that depend on mi and TO2. In a spatially varying mi- 
crowave field, this results in a state-dependent potential 
landscape. In addition, there is a small common mode 
part of the potential due to the electric field of the mi- 
crowave, which shifts all levels \F, mp) identically. 

A trap for neutral atoms based on microwave po- 
tentials has been proposed in and experimentally 
demonstrated in [l!j. This trap employs microwave ra- 
diation in the far field of the source. Unlike the case of 
optical radiation, which can be tightly focussed due to 
its short wavelength, the centimeter wavelength A mw of 
microwave radiation poses severe limitations on far-field 
traps: field gradients are very weak [l9| and structuring 
the potential on the micrometer scale is impossible. 

On atom chips, there is a natural solution to this prob- 
lem. The atoms are trapped at distances d -C A mw from 
the chip surface. Thus, they can be manipulated with 
microwave near fields, generated by microwave signals 
in on-chip transmission lines |20| . In the near field of 
the source currents and voltages, the microwave fields 
have the same position dependence as the static fields 
created by equivalent stationary sources. The maximum 
field gradients depend on the size of the transmission 
line conductors and on the distance d, not on A mw . In 
this way, state-dependent microwave potentials varying 
on the micrometer scale can be realized. In combina- 
tion with state-independent static magnetic microtraps, 
the complex potential geometries required for QIP can 
be realized. 

The use of microwaves to create the state-selective po- 
tential is essential for our proposal, since a combination 
of static magnetic and electric fields, as considered in 
[ToL l2l| . does not provide state-selective potentials for 
our robust qubit state pair, whose magnetic moments and 
electrostatic polarizabilities are equal. Optical potentials 
created by focussed laser beams with a frequency close 
to the Dl or D2 transition of 87 Rb are also impractical: 
if the detuning of the laser from the atomic resonance 
is much larger than the hyperfine splitting of the 87 Rb 
ground state, the resulting optical potentials are again 
nearly identical for the states |0) and |1). If, on the other 
hand, a detuning comparable to the hyperfine splitting 
is used, a differential optical potential could be created, 
but problems with decoherence due to spontaneous scat- 
tering of photons arise. 
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m F = -2 m F = -1 m F = m F = 1 m F = 2 



FIG. 2: (Color online) Energy-level diagram of the hyperfine 
structure of the 87 Rb ground state in a combined static mag- 
netic and microwave field. y^' mp indicates the energy shift 
due to the static Zeeman effect, which is (nearly) identical 
for |0) and |1). Due to the magnetic field of the microwave, 
the level \F,m F ) is shifted in energy by li*" F . This shift 
has opposite sign for |0) and |1), as indicated in the figure 
(the shift of the other levels is not shown). The microwave 
transitions contributing to the shift of |0) and |1) are shown 
for A^'™ 2 > (blue detuning). 



A. Theory of microwave potentials 

We consider the hyperfine levels \F,mp) of the 551/2 
ground state of a 87 Rb atom in a static magnetic field 
Bo(r). In addition, the atom is exposed to a microwave 
of frequency lu with magnetic field B mw (r) cos(wi) and 
electric field E mw (r) cos(wt + ip), where ip accounts for 
a phase shift between electric and magnetic field. We 
now calculate the adiabatic trapping potentials for the 
atomic hyperfine levels in the combined static and mi- 
crowave fields. The trapping frequencies considered in 
this paper are sufficiently small so that non-adiabatic 
transitions due to the atomic motion do not limit the 
trap lifetime 0. 

The static magnetic field Bo(r) leads to a Zeeman shift 

Vz' mF (r) = hb9f m F |B (r) (3) 

of the levels (see Fig.[2J), where /is is the Bohr magneton 
and gF is the hyperfine Lande factor for the state F. This 
shift is identical for states |0) and |1) and we identify it 
with the common potential in (J2J, 

u c (r) = F|< +1 (r) = Vh>-\r) = ^|B (r)|. (4) 

Effects of the nuclear magnetic moment have been ne- 
glected in |(2J), they lead to corrections in the potential of 
the order of 10 -3 . Here and in the following, we take the 
local direction of the static field Bo(r) as the quantiza- 
tion axis for the levels \F,mp), 



The magnetic component of the microwave field cou- 
ples the |l,mi) to the \2,m2) sublevels, with Rabi fre- 
quencies 

2,m 2/ ^_ (2,TO 2 |A-B mw (r)|l,TOi) 

for the different transitions. We apply the rotating-wave 
approximation, which introduces negligible error for the 
parameters considered in this paper. In (JSJ, fx = /isgjJ 
is the operator of the electron magnetic moment (gj ~ 2); 
its matrix elements are evaluated in the basis \F,mp) 
defined with respect to the local quantization axis along 
Brj(r). In a combined static magnetic and microwave 
trap, as considered here, both B (r) and B mw (r) vary 
with position. This leads to a position-dependent mi- 
crowave coupling with in general all polarization compo- 
nents present. The detuning of the microwave from the 
resonance of the transition |I, mi) — > |2, 7712) is given by 

A^(r) = A -^|(m 2 +m 1 )|Bo(r)|, (6) 

where Ao = oj — ujq is the detuning from the transition 
|F = 1,to f = 0) -> \F = 2,m F = 0), and the dif- 
ferent Zeeman shifts © of the levels have been taken 
into account. In the following, we concentrate on the 
limit of large detuning |A^| 2 > |«i;™J| 2 , which al- 
lows for long coherence lifetimes of the qubit states in 
the microwave potential (see Sec. IVIlT^I . In this limit, 
the energy shifts due to the microwave coupling can be 
evaluated perturbatively for each transition. The over- 
all magnetic microwave potential for the level \F,rriF) 
equals the sum of the energy shifts due to the individual 
transitions connecting to this level. Thus, the magnetic 
microwave potentials for the sublevels of F = 1 are given 
by 

h Ifi 2 '™ 2 ^! 2 
y i,m 1(r] = ^ylWi (7) 

"mw \ L ) , A 2,m 2 / s ' V > 

while the potentials for F — 2 are given by 

K 2 *( r ) = --Vl^™i^ ( 8 ) 

mi a l,miW 

As desired, the potentials for F = 1 and F — 2 have 
opposite sign, leading to a differential potential for the 
qubit states |0) and In Fig. [21 the relevant transi- 
tions contributing to the potentials for |0) and |1) are 
shown for a general microwave field with all polarization 
components present. From Fig. it is evident that the 
microwave potentials for the two states will not have ex- 
actly the same position dependence (in addition to the 
difference in sign), since only two polarization compo- 
nents contribute to the potential for |1), while all three 
components contribute to the potential for |0). Further- 
more, by reversing the sign of the detuning A^'™J for all 
levels, the overall sign of the microwave potentials and 
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therefore the role of the states during the gate operation 
can be interchanged, i.e. either |0) or |1) can be chosen 
as the oscillating state. We choose A 2 '™^ > 0, resulting 
in |1) as the oscillating state, since the potential for |1) 
is closer to harmonic in our actual trap design. 

In addition to the magnetic microwave field, the elec- 
tric field E mw (r) cos(cut + cp) also leads to energy shifts. 
The electric field of the microwave polarizes the atoms, 
leading to a time-averaged quadratic Stark shift 

K 1 (r) = -||E mw (r)| 2 , (9) 

where a is the scalar DC polarizability of the 87 Rb 
ground state. In 0, we have averaged over the fast os- 
cillation of the microwave at frequency lo, which is much 
faster than the atomic motion. Eq. J5J is actually the 
DC limit of optical dipole potentials . Since all levels 
\F,mp) belong to the electronic ground state, their po- 
larizabilities are equal, leading to an electric microwave 
potential which is the same for all levels. 

The total microwave potential for state |0) is identified 
with uo(r) in @, 

uo(r) - V el (r)+V£-\r) (10) 

_ « 2 ft ° 

while the microwave potential for state |1) is identified 
with ui(r) in (J2J, 

Ul (r) = V el (r) + V£+\r) (11) 
- -4 |Emw(r)l ~ 4 ^ A 2 ' +1 (r) ' 

mi =0 ^l,mi V r J 




FIG. 3: (Color online) Chip layout for the microwave colli- 
sional phase gate, (a) Cut through substrate, (b) Top view 
of wire layout. Wire parameters: w = 0.8 pm, s = 0.1 pm, 
t — 0.2 (j,m, dr = 4.2 /im, u>t = 1.5 fim, tr ~ 1.0/im. The 
lower wire (current It) is fabricated into a groove which was 
etched into the substrate, and covered by an insulating layer. 
The DC currents It, Ic, II, and In and the orientation of 
the magnetic bias field B used to create the state-independent 
potential u c (r) are shown. The atomic wave functions in this 
potential are indicated. The three upper wires form a mi- 
crowave coplanar waveguide (CPW), compare Fig. 0] 



IV. MICROWAVE CHIP DESIGN 

To realize the state-selective potential for the atoms, 
we propose to use a microwave atom chip with a wire lay- 
out as shown in Fig. El The chip consists of two layers of 
gold metallization on a high resistivity silicon substrate, 
separated by a thin dielectric insulation layer. The wires 
carry stationary (DC) currents, which, when combined 
with appropriate stationary and homogeneous magnetic 
bias fields, create the state-independent potential u c (r). 
In addition to carrying DC currents, the three wires on 
the upper gold layer form a coplanar waveguide (CPW) 
[2fH | for a microwave at frequency ui. The microwave fields 
guided by these conductors create the state-dependent 
potential Ui(r). The combination of DC and microwave 
currents in the same wires is possible by the use of bias 
injection circuits [23l |. 

A. Static magnetic trap 

For the phase gate with three oscillations (N = 3, see 
Sec. IVI|) . an initial static trap with DC currents It — 



348.440mA, I c = -0.813 mA, and I L = I R = 1.204 mA 
is used (see Fig. |3J). The components of the homoge- 
neous bias field B = (B x ,B y ,B z ) are B x = -4.464G, 
B y = 103.717 G, and B z = 0.000 G. For assumptions 
about the stability and accuracy of the currents and mag- 
netic fields, see Sec. IVII Fl In the calculation of the po- 
tential, the finite size of the wires is taken into account. 
For the parameters given above, u c (r) is a Ioffe-type dou- 
ble well potential |2| along an axis x' in the xy-plane, 
which is tilted by a small angle 9 = 0.02 with respect to 
x. The distance of the double well from the wire sur- 
face is d — 1.80 /im. The magnetic field in the trap 
center is |E$o| = 3.230 G, which maximizes the coher- 
ence time in the absence of microwave coupling [24| . B 
is directed approximately along —x. The distance be- 
tween the minima of the double well is d x — 1.32 /im. In 
the transverse dimensions, the trap provides a tight har- 
monic confinement with almost identical trap frequencies 
ujy/2ir — oj z /2it = oj±/2it = 77.46kHz. For the simula- 
tion, we determine uj±(x) as a function of the longitudinal 
coordinate x, although the relative variation along x is 
only < 1 x 10~ 3 . The axial trap frequency of the poten- 
tial wells is ui x /2ir = 4.432 kHz. All traps considered in 
our gate simulation satisfy uj x <C oj±_ ■ 
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B. Microwave propagation on the chip 

To accurately determine the state-selective part of the 
potential, a simulation of the microwave electromagnetic 
fields in the proximity of the conductors transmitting the 
microwave signal is necessary. Since the transverse {xz- 
plane) size of the CPW is only a few micrometers, much 
smaller than the centimeter wavelength of the microwave, 
a quasi-static field analysis is sufficient to determine the 
microwave propagation characteristics |20| . At oj/2tt = 
6.8 GHz, the skin depth is S s = ^/2/uifio<r — 0.9 /im, 
where a = 4.5 x 10 7 ft _1 m _1 is the conductivity of gold. 
The skin depth has to be compared to the thickness t 
of the CPW conductors, see Fig- El Since t < 6 S , the 
microwave electromagnetic field is not screened from the 
inside of the conductors, and microwave currents flow in 
the whole cross section of the gold wires, not only in 
a thin layer below the conductor surface. Taking this 
into account is equivalent to a proper treatment of con- 
ductor losses in the microwave simulation. The common 
assumption of perfect conductor boundary conditions on 
the wires is no longer justified. The quasistatic charac- 
teristics of a micron-sized CPW with a proper treatment 
of conductor-loss effects are given in [25|, |2(j • Compared 
to conductor losses, dielectric losses are negligible if a 
high-resistivity Si substrate is used, which is a standard 
substrate material for microwave circuits. 

We have performed a quasistatic analysis of the CPW 
in Fig. neglecting the wire in the lower plane of met- 
allization. The analysis consists of an electrostatic and 
a magneto-quasistatic simulation of the fields in a two- 
dimensional transverse (xz-plane) cross section of the 
waveguide. The results of this simulation are shown 
in Fig. 0] From the fields, the complex characteristic 
impedance Z c of the CPW can be determined. For our 
CPW, Z c = 130 ft ■ e - l - 247r . The microwave propaga- 
tion constant /3 mw = 27r/A mw and the attenuation con- 
stant a mw are of comparable magnitude, /3 mw ~ a mw ~ 
1 x 10 3 m~ 1 at 6.8 GHz. Damping and propagation ef- 
fects are negligible on a micron scale. The influence of the 
lower wire on the microwave propagation characteristics 
of the CPW is expected to be negligible, since dr ^> s. 
To estimate the effect of the lower wire, we have included 
a homogeneous gold layer in the plane of the lower wire 
in our two-dimensional simulation. This leads to a small 
relative change in |Z C | of 1 x 10~ 2 . 

The microwave signal on the CPW is given by a com- 
plex voltage amplitude V on the center conductor, while 
the two outer conductors are set to V = (ground), see 
Fig. 0] Here we use phasor notation to describe the mi- 
crowave signal and suppress the harmonic time depen- 
dence e luJt (2£j. The amplitude of the microwave cur- 
rent in the center conductor is Io, while the two outer 
conductors carry a current — 1q/2. Voltage and current 
are related by Vo = Z c Io, i.e. there is a phase shift 
(p = arg(Z c ) between them, which gives rise to the phase 
shift <p between electric and magnetic field of the mi- 
crowave. The microwave potentials, however, are inde- 



pendent of (p (see Sec. IIII A"|) . 

In the simulation of the gate dynamics and in Fig. 
the amplitudes of the microwave signal on the coplanar 
waveguide are 

| Vol = y/Xtfj • 1.9&95V (12) 
|I | = iJ4 = V^W- 1 5-343mA (13) 

l Z c| 

Here, A(i) is the modulating function of the microwave 
potential during the gate operation, which will be op- 
timized by optimal control techniques as described in 
Sec. IVII The detuning of the microwave from the tran- 
sition \F = l,m F = 0) -> \F = 2,m F = 0) is 
Aq = 2tt ■ 29.4 MHz. For these parameters, the ratio 
between the electric and the magnetic microwave poten- 
tial is Vci ~ 0.2Vi nw , i.e. the state-selective magnetic 
part of the microwave potential is dominating. This is 
important since V^i reduces the barrier of the double well 
also for state |0). It is possible to avoid this unwanted 
effect of V e \ by a small adjustment of the magnetostatic 
potential during the gate operation. To compensate V e i, 
we modulate 

B x {t) = (—4.464 + Ao(i) ■ 0.036) G (14) 
I c (t) = (-0.813- A (t) -0.039) mA, (15) 

where Aq(£) ~ X(t). The function Ao(i) consists of simple 
linear ramps, as explained below fSec lVTA")! . 

If the microwave is turned on to full power (A — 1 and 
Ao = 1), the barrier of the static double well trap is re- 
moved for state |1), leaving a single potential well of lon- 
gitudinal trap frequency W| 1 ^/27r = 5.448 kHz, as shown 
in Fig. 2] For state |0), the two potential wells are shifted 
apart and the trap frequency in each of the wells changes 
to cj| )/27T = 4.775 kHz. The transverse confinement of 
the static trap is unchanged, both d and lo \_ change by 
less than 10 -3 . Although we include the full position 
dependence of itj(r) in the simulation, it would be suffi- 
cient to consider the state-selective potential Ui(x, yo, zq) 
at the transverse position of the static trap minimum 

(yo,z ). 

C. Chip fabrication 

The fabrication of the chip relies on electron beam 
lithography, which provides sub-micron resolution, in 
combination with lift-off metallization techniques and 
anisotropic etching of the Si substrate. The lower gold 
wire is fabricated as follows: prior to the metal deposi- 
tion, a groove is etched into a high-resistivity Si substrate 
by anisotropic reactive ion etching, using lithographically 
patterned photoresist as etch mask. The gold wire is 
subsequently fabricated into this groove using the same 
photoresist for lift-off. The remaining groove is filled and 
the wire is insulated by depositing a dielectric layer, such 
as SiC>2 or polyimide. Since the lower wire is contained 
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x \pm] x [fim] 

FIG. 4: (Color online) Transverse (:rz-plane) magnetic and electric fields of the microwave on the coplanar waveguide. The 
fields were obtained by a quasi-static simulation including conduct or- loss effects. The longitudinal (y-direction) electric field, 
which is several orders of magnitude smaller than the transverse electric field, is not shown in the plots, (a) Magnitude of the 
fields for A = 1 shown in a cross section of the CPW. (b) Transverse microwave field components as a function of a; at a distance 
d — 1.80 /im from the wire surface, corresponding to the line z — zo — 1.90 /im in (a), which indicates the z-position of the 
static trap minimum, (c) Microwave voltage and current amplitudes on the CPW used in the simulation, a is the conductivity 
of the gold wires, e r the dielectric constant of the Si substrate. For the size of the conductors, see Fig. EH 



in the groove, this insulating layer can be very thin and 
still provide good planarization of the surface. In this 
way, good thermal contact of the wires to the substrate 
is guaranteed. On top of the insulating layer, the CPW is 
fabricated using standard electron-beam lithography and 
lift-off. 



V. SIMULATION OF THE GATE OPERATION 

The gate operation is simulated numerically by solv- 
ing a time-dependent Schrodinger equation of the two- 
particle dynamics along the direction of the double well 
x 1 (in the following, we write x instead of x 1 to sim- 
plify notation) . We assume that the atoms remain in the 
transverse ground state of the trap throughout the gate 
operation; this assumption is justified in Sec. IVIBI and 
Sec. IVII Bl The initial state of the two atomic qubits in 
an arbitrary superposition of internal states \ij) is 



|*(t = 0)> = 



v*,je{o,i} 



ot.ij\ij) ® \ip(xi,x 2 ,t = 0)), 



(16) 

where \ij)(xi,x 2 ,t = 0)) describes the initial motional 
state of the two atoms, which is independent of the in- 
ternal state, with one atom in each well of the potential 
u c (x) (see Fig.^i). During the gate operation, the inter- 
nal and motional states of the atoms are entangled due 



to the dynamics in the state-selective potential: 

!*(*)) = "tfly)® \^a(xi,x 2 ,t)}. (17) 

i,j'e{o,i} 

The dynamics of \ipij(xi,x 2 ,t)) (see Fig.^i) is governed 
by the time-dependent Hamiltonian 

Hij(t) = T Xl +T X2 +Ui(xi,t)+Uj(x2,t)+Vy t (\x 2 -Xi\,t), 

(18) 

where T Xi denotes the kinetic energy operators and 



V&{\X2-Xt\,t) 



1 - lAQal 3 /a±(x 1 ,t) 



S(\x 2 - Xl \) (19) 



is the effective one-dimensional interaction potential be- 
tween the atoms a\ J ~ 5.4 nm is the s-wave scatter- 
ing length for collisions of 87 Rb atoms in state \ij), and 
a±(x%,t) = y/2h/muj±(xi, t) is the size of the ground 
state in the transverse direction. The transverse trap fre- 
quency u>j_(xx,t) can be used as a control parameter to 
optimize the dynamics and is therefore time-dependent. 
The numerical simulation of the dynamics is performed 
using a Fast Fourier Transform method. 

The entanglement between the motional and internal 
states of the atoms is crucial for establishing the col- 
lisional phase shift in state |11). However, at the end 
of the gate operation, the motional and internal states 
have to factorize again, so that the truth table (JTJ is 
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implemented. Any residual entanglement between the 
motional and internal states leads to a reduction of the 
gate fidelity. To avoid this, the wave function of the 
atoms has to show a revival at t — r g and regain its ini- 
tial form, apart from the phase factor in state |11). To 
quantify this, we define the overlaps 

Oij(t) = (i) ij (x 1 ,x 2 ,t)\i;(x 1 ,X2,Q)), (20) 

from which the overlap fidelities 

F<i® = \Oij(t)\ 2 (21) 

are calculated. The evolution of the gate phase is given 
by 

4> g (t) = (f>n(t) + 0oo (*) - 0oi (t) - 0io (t) , (22) 

where we have defined the collisional phases 

4>ij(t) = axg[(ipij(xi,X2,t)\ili^(xi,X2,t)}] . (23) 

Here, \ipfj(xi,X2,t)) is the motional state evolved from 
\ip(xi, X2, 0)) in the absence of collisions ^(j- fyijiP) is 
well defined even for times t where Oij(t) ~ 0. Optimal 
gate performance corresponds to Fij(r g ) = ljVjj, i.e. a 
complete revival of the wave function, and a gate phase 
4>g( T g) = 7r - In a realistic non-harmonic potential, how- 
ever, the revival of the atomic wave function will be in- 
complete, and the collisions between the atoms will lead 
to an additional distortion of the wave function. Both ef- 
fects will lead to a reduction of the gate fidelity, and have 
to be avoided as much as possible by an optimization of 
the gate dynamics. 

VI. OPTIMIZATION OF THE GATE 

We apply optimal control techniques (2^, |3(j to op- 
timize the gate performance. This is performed in two 
steps. First, the dynamics without atom-atom interac- 
tions is optimized using X(t) as control parameter. In this 
way, we optimize the revival of the wave function in the 
non-harmonic potential. In a second step, the dynamics 
in the presence of interactions is optimized further using 
oj±(t) as control parameter. This reduces wave packet 
distortion due to collisions. The two steps are discussed 
in the following subsections for a gate with N = 3 oscil- 
lations. 



A. Optimal control of A(i) 

In the first stage of optimization, X(t) is used as control 
parameter. Experimentally, this corresponds to control- 
ling the microwave power. We start with a trial function 
Ao(i) which consists of linear ramps, as shown in Fig. [SJu. 
This choice represents a compromise between a sudden 
removal of the barrier for state |1) and an adiabatic shift 



of the potential wells for state |0). We optimize r g and 
the slope of the linear ramps of Ao (t) by hand to achieve 
an initial gate fidelity > 0.95. Starting with A(t) = Ao(t), 
the function X(t) is then adjusted by an optimal control 
algorithm, which neglects atom-atom interactions in or- 
der to provide faster convergence. This corresponds to 
setting a s = in l|19|) . In the absence of interactions, 
it is sufficient to optimize the single-particle overlaps 
Oi = \{ijji(x,Tg)\ijj(x,0))\, which are optimized simulta- 
neously for i € {0, 1} by the algorithm using the Krotov 
method After the optimization has converged, the 
resulting optimized X(t) is used in a simulation of the two- 
particle dynamics with interactions (a s ^ 0), in order to 
determine the two-particle overlap fidelities Fij(r g ). The 
strength of the interaction (|19fl was adjusted by setting 
the transverse trap frequency u>± to a value constant in 
time such that the gate phase is <j) g ~ n. The result of 
this first stage of optimal control is the microwave power 
modulation function X(t) shown in Fig. [SJi, which devi- 
ates from Ao(t) by small modulations, and the fidelities 
shown in the left part of Table Q] 

The trial function Xo(t) is also used to modulate the 
magnetostatic potential in order to compensate the effect 
of the microwave electric field, as explained in Sec. IIVBI 
We choose Ao(t) instead of X(t) for this modulation in 
order to keep the number of experimental parameters 
which are subject to optimal control as small as possible. 

Our simulation shows that collisional interactions be- 
tween the atoms are negligible in basis states other 
than 1 11), with 0oo(t 3 ) = and 0oi(t 5 ) = 0io(t 3 ) ~ 
10~ 3 0n(r g ). In Fig. [Sji, the phase evolution during the 
gate operation is shown. The main effect of the inter- 
action is to provide the phase shift in state |11), corre- 
sponding to the steps in 4> g {t) each time a collision in 
this state takes place. An undesired effect of the colli- 
sion is to decrease the fidelity fn(r 9 ) compared to the 
non-interacting case due to wave packet distortion dur- 
ing the collision. We have performed the optimization 
of the gate for different numbers of oscillations N. For 
each value of N, a different value of uj± was chosen to 
adjust the gate phase to <j) g ~ tt. For smaller N, the 
gate phase has to be acquired in a smaller number of os- 
cillations, therefore the interactions need to be stronger 
and the phase shift per collision is larger. Correspond- 
ingly, the collisional distortion of the wave function is also 
larger, and Fn is smaller. This general tendency can be 
seen in the left part of Table [I] by comparing the values 
of Fix for different N. The best gate performance was 
achieved for N = 5, with Fn = 0.991 and a gate time of 
r g = 1.838 ms. 



B. Optimal control of wj_(t) 

After optimization of A(i), the fidelity of the gate is 
limited by collisional distortion of the wave function in 
state |11). To overcome this limitation, we have imple- 
mented a second stage of optimization, in which the inter- 
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FIG. 5: (Color online) Dynamics during the gate operation, shown for iV = 3 oscillations, the gate time is r 9 = 1.110 ms. 
Both X(t) and LO±(t) are used as control parameters, (a) Optimal control of the microwave power. Upper plot: Initial trial 
function Xo(t) (dashed) and optimized control parameter X(t) (solid line). Each triangular ramp of X(t) corresponds to a full 
oscillation of state |1). Lower plot: the difference SX(t) — X(t) — Xo(t) shows small modulations, (b) Optimal control of the 
effective one-dimensional interaction strength via modulation of the transverse trap frequency uj±(t). (c) Evolution of the 
overlap fidelities during the gate operation: Foo(t) (dashed), Foi(t) = Fio(i) (dotted), Fn(i) (solid line), (d) Evolution of the 
gate phase <j> g (t). The phase shift steps are due to the six collisions in state |11). 



Optimal control of X(t) 



Optimal control of X(t) and ui±(t) 



N r g [ms] 


Oo Oi Foo foi Fn (j>g j 7r 


Foo F i Fn <f>g/iv 


2 0.696 

3 1.110 

4 1.389 

5 1.838 

6 2.219 


0.996 0.996 0.993 0.993 0.960 0.996 
0.998 0.997 0.996 0.996 0.986 0.999 
0.996 0.995 0.991 0.991 0.991 0.991 
0.999 0.997 0.997 0.997 0.991 0.995 
0.996 0.998 0.992 0.993 0.985 0.993 


0.993 0.993 0.987 0.998 
0.996 0.996 0.995 0.997 



TABLE I: Optimized gate performance for different numbers of oscillations N and correspondingly different gate times r fl . All 
two-particle fidelities Fij, single-particle overlaps Oi, and the gate phase (j> g are evaluated at t = T g . Either only X(t) (left part 
of the table) or both X(t) and u±(t) (right part of the table) were chosen as control parameters. The optimization of u>±(t) was 
performed only for N < 3, which is particularly interesting due to the short gate times r 9 < 10~ 3 r t . 



action potential (|19fl is controlled during the gate opera- 
tion with the transverse trap frequency uj± (t) as control 
parameter. Experimentally, this corresponds to control- 
ling the bias field component B y (t) and the DC current 
2t(£) proportional to the desired modulation of ui±(t). 
Since the collisions in internal states other than |11) 
are negligible, this optimization affects only state |11). 
We optimize the overlap fidelity Fn(r s ) and the phase 
4>ii(r g ) simultaneously, using again the Krotov method. 
For the microwave power we use the optimized X(t) de- 
termined in the previous section. 

Care has to be taken that the modulation of uj±(t) 
does not create excitations of the transverse state of the 
atoms. Since it is not possible to put constraints on the 
time derivative of the control parameter in the Krotov 
method of optimization, we have chosen the following 
strategy to avoid transverse excitations. We parameter- 



ize the transverse trap frequency 

u±(t) = w_L(0)[Atanha(t) + 1], (24) 

with a(t) being the dimcnsionlcss optimal control pa- 
rameter, and Cl>j_(0) is determined for a given N as in 
the previous section. In this way, the optimization is 
constrained to a maximum modulation amplitude set 
by A. After the optimization on a(t) has converged, 
we remove high-frequency components from the result- 
ing modulation of uj±(t) by filtering with a cutoff fre- 
quency w c . We choose lu c < 2luj_(0) in order to avoid 
parametric excitation of the transverse degrees of free- 
dom. By performing the gate simulation again with 
the filtered modulation LUj_(t), we determine the over- 
lap fidelities and the gate phase. We have seen that 
the filtering does not significantly decrease the fidelity 
if lo c 3> uj x . This shows that the high frequency com- 
ponents with to > lo c were artifacts of the optimization 
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algorithm without physical significance. The modulation 
of u>±(t) after filtering is shown in Fig.03. For this mod- 
ulation, A — 0.2 and oj c = 0.8u;_i_(0). To quantify the 
transverse excitation probability p™(i) during the gate 
operation, we simulate the transverse dynamics using a 
harmonic oscillator model, where the frequency is mod- 
ulated according to the filtered u>±(t). This simulation 
yields p™(t) < 7 x 10 -4 , therefore transverse excitations 
due to the modulation of u>±(t) do not limit the gate 
performance. 

Fig. [S] is the main result of our simulation. It shows 
the gate dynamics for N = 3 after optimization of both 
X(t) and u±(t). Corresponding numbers for the overlap 
fidelities F%j and the gate phase <p g are shown in the right 
part of Table [fl The improvement compared to optimal 
control of X(t) alone (left part of Table P) is an increase 
of F n from F xl = 0.986 to F n = 0.995, which is now 
comparable to the overlap fidelities in the other states. 
The gate time is r g = 1.110 ms. 

C. Gate fidelity 

In order to provide an estimate of the overall gate per- 
formance, we use the definition given in [lCj for the gate 
fidelity: 

F = min x {Tr cxt [( X \U S [\ X )(x\ ® Po) SW\ X )] } ■ (25) 

Here, |x) is an arbitrary internal state of both atoms, and 
\x) is the state resulting from \x) using the actual trans- 
formation. U and S are the operators for time evolution 
and symmetrization under particle interchange, respec- 
tively. The density matrix po is the initial two-particle 
motional state. According to this definition, the gate 
fidelity is 

F = 0.997 for N = 3, r g = 1.110ms (26) 

after optimization of both X(t) and u>±(t). Even after 
the optimization, wave packet distortion still contributes 
the largest error reducing the fidelity. The error sources 
discussed in the next section lead to an additional error 
of the order of 1 x 10~ 3 , mainly due to the finite trap 
lifetime, which is limited by surface effects. If we include 
this error in the calculation of the fidelity, we get F — 
0.996. By comparison, with trap and coherence lifetimes 
Tt ~ t c ~ 10 3 T g , the maximum achievable fidelity of the 
gate with N = 3 would be F = 0.999, if wave packet 
distortion due to interaction dynamics could be reduced 
to negligible values. 

VII. ERROR SOURCES 

In this section, we discuss several effects which could 
possibly limit the fidelity of the phase gate proposed here, 
and justify the assumptions made in the description of 
the gate dynamics. 




k B T/huj x 

FIG. 6: Gate infidelity 1 — F(T) as a function of temperature 
for the gate with N = 3. 

A. Finite temperature 

The result for the fidelity i)26|) is obtained for zero tem- 
perature (T = 0) of the atoms. This corresponds to an 
initial motional state with one atom in the ground state 
| ni = 0} of the left potential well and the other atom 
in the ground state \n 2 = 0) of the right potential well 
of u c (r) . To study the effect of finite temperature T on 
the gate performance, we start with an initial density 
operator for the motional degrees of freedom 

oo 

Pa = p nx,nAT)\ni){ni\®\n 2 ){n 2 \, (27) 

where P ni)Tl3 (T) is the probability for the occupation of 
the motional state \rii) for atom 1 and \n 2 ) for atom 2 in 
the initial double well trap. The probabilities P ni ,n 2 (T) 
are calculated assuming a thermal distribution corre- 
sponding to a temperature T in the canonical ensemble. 
We have evaluated the fidelity F = F(T) of our gate 
with N = 3 as a function of T, see Fig. HO The results 
show that high fidelity gate operations are only possible 
if the two atoms can be prepared in the motional ground 
state of the trap, with very low occupation probability 
of excited states. We assume that kBT/huj x < 0.1 can 
be reached, corresponding to a temperature T < 20 nK 
in the initial double well trap. In this case, the fidelity is 
not reduced significantly due to the finite temperature. 

B. Condition for quasi-lD dynamics 

To ensure one-dimensional dynamics of the atoms in 
the experiment and to justify the one-dimensional simu- 
lation of the gate presented above, transverse excitations 
due to the collisions in state |11) have to be avoided. This 
mechanism could limit the fidelity even if transverse ex- 
citations due to the modulation of the transverse trap 
frequency co±(t) are suppressed (cf. Sec. I VI Bl . 

We calculate the probability (t) of transverse excita- 
tions due to collisions using a simple model. We consider 
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FIG. 7: Transverse excitation probability p c ±(r g ) due to col- 
lisions of the atoms for N = 3. p±(T g ) is shown as a function 
of the kinetic energy Skin of one atom in state |1) at the time 
of the collision. 



two interacting atoms in a three-dimensional anisotropic 
harmonic oscillator potential, with u> x <C uj v = lu z = luj_. 
Exact solutions for this model are known |22j ■ The atoms 
are initially in the transverse ground state of the trap, 
but separated in two harmonic potential wells with sep- 
aration d x along the x direction, similar to the initial 
state of our phase gate. At time t = 0, the initial double 
well is switched off and the atoms evolve in the single 
anisotropic harmonic oscillator potential. The dynamics 
is obtained by expanding the initial three-dimensional 
two-atom state on the interacting basis given in |32( . For 
each evolution time t, we compute the reduced density 
matrix of the transverse motion by tracing out the axial 
degrees of freedom. In this way we compute the prob- 
ability p c j_(t) as a function of time. We have performed 
this calculation for different d x , corresponding to differ- 
ent kinetic energies -Ekin of one atom at the time of the 
collision. In Fig. 0we show Px( T s) f° r N — 3 as a func- 
tion of 2 Eki n /hu>±. We find that transverse excitations 
are energetically suppressed if 



2 E kia < Tiuij 



(28) 



For all trapping geometries used in our simulation of the 
phase gate, 2 E^/Tllu^ < 0.7, and the transverse excita- 
tion probability is suppressed to negligible values. 



C. Microwave coupling and qubit dephasing 

In contrast to optical dipole transitions, spontaneous 
emission is negligible for microwave transitions between 
ground state sublevels. Therefore microwaves can also be 
used to create adiabatic potentials with the microwave 
frequency tuned to resonance with a particular hyper- 
fine transition, similar to the adiabatic radio-frequency 
potentials investigated in [!, US 13 • However, resonant 
coupling results in strong mixing of the hyperfine levels 
connected by the microwave. This is undesirable for our 
qubit state pair, since it would destroy the excellent co- 
herence properties of the qubit by admixtures of other 



states with different magnetic moments. For this reason, 
we concentrate on the limit of large microwave detun- 
ing in this paper. The large detuning also ensures that 
the modulation of the microwave power during the gate 
operation is adiabatic with respect to the internal state 
dynamics of the atoms. 

For the trap parameters considered here, 



max — < 
A^ 



itr 



(29) 



where we have set flu = fl{ m \ and A = A 1 ' m ^ to 
shorten notation in this section. In the absence of the 
microwave coupling, the differential magnetic moment 
Sfi = dE\\)/dB — dE\ }/dB of the qubit states |1) and 
|0) can be calculated from the Breit-Rabi formula [24| . 
For magnetic fields Bq < 6 G, we get \5fi\ < 2 x 10~ 3 /is, 
with 5/i = 0, i.e. a vanishing first-order differential Zee- 
man shift, for Bq = 3.229 G |24|, as in the center of our 
static trap. In the presence of the microwave coupling, 
5fi changes due to the admixture of other magnetic sub- 
levels, whose magnetic moments differ by multiples of 
^b/2 from the magnetic moment of the qubit states. For 
ft 2 R <C A 2 , the order of magnitude of this change can be 
estimated to 



5jJL 



O 2 



4A 



2 Hb < 2.5 x 10 3 /j, B , 



(30) 



where we have used (|29|) . 

The differential magnetic moment determines the co- 
herence time t c of the qubit in the presence of longitu- 
dinal magnetic field fluctuations (pure 'dephasing' of the 
qubit). The frequency spectrum of the fluctuations plays 
an important role in this context: Low frequency fluc- 
tuations of the mag netic field are most harmful to the 
coherence time j35(, since they do not average out on 
the time scale T g of a single gate operation. In contrast 
to magnetic near-field noise arising from the chip surface 
(see Sec. lVIIDj) . which has a flat spectrum in the relevant 
frequency range |15| . technical magnetic field noise typi- 
cally increases towards low frequencies and dominates in 
the frequency range < u < 10/r g considered here. 

The atom chip is surrounded by magnetic shielding, so 
that the residual fluctuations of the magnetic field are 
limited by the stability of the current sources used for 
the wires and bias magnetic field coils. We assume that 
fluctuations in the frequency range < u> < 10/r g can 
be reduced to an r.m.s. amplitude of SB < 0.01 mG, 
consistent with the stability of the current sources in our 
experiments. Under this assumption, we expect a coher- 
ence time of the order of 



5fi SB 



5 s, 



(31) 



where we have used (13011 . Comparing this value of t c with 
the trap and coherence lifetimes due to surface effects 
calculated in Sec. IVH Dl we find that pure dephasing of 
the qubit does not limit the fidelity of the gate. 
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D. Trap loss and decoherence due to the chip 
surface 



A fundamental source of trap loss, heating, spin- and 
motional decoherence on atom chips is magnetic near- 
field noise originating from thermal currents in the chip 
wires |Tfjj . For our state pair, pure spin dephasing arising 
from longitudinal magnetic near-field noise is negligible. 
The dominating surface effect is trap loss due to spin-flips 
induced by transverse magnetic field noise. For the chip 
layout shown in Fig. [3] and a distance d = 1.80 /im from 
the chip surface, we estimate an average spin-flip rate of 
r s = 0.9 s _1 , taking the finite thickness and width of the 
wires |3(| and the different matrix elements for spin flips 
of state |0) and |1) into account. This corresponds to a 
trap lifetime of 



n = r- 1 = lis. 



(32) 



Compared with surface effects, loss rates due to collisions 
with background gas atoms are negligible. The rates for 
heating and motional decoherence due to magnetic near- 
field noise are comparable to r s ^J, so that we expect 
an overall coherence lifetime r c ~ r t . For N = 3, the gate 
operation time is r g — 1.110 ms. The mentioned surface 
effects introduce an error of 1 — exp(— T s r g ) = 1 x 10~ 3 in 
the gate operation. This error is smaller than the error 
due to wave packet distortion, and therefore does not 
significantly decrease the gate fidelity (|26l) . The error 
could be further reduced by reducing the wire thickness 
or using a trap at a slightly higher atom-surface distance. 



the wires is approximately constant over the whole wire 
cross section. Therefore the resistance of the wires is 
given by the DC resistance and the maximum tolerable 
microwave current density is comparable to the maxi- 
mum DC current density in wires of this size. In our 
simulation, we assume total (microwave + DC) current 
densities jtot < 1 x 10 11 Am -2 in the wires forming the 
coplanar strip line. The DC current density in the lower 
wire is jtot = 2 x 10 11 Am" 2 . Comparable current den- 
sities have been realized experimentally in |8j] on Si sub- 
strates with a 20 nm Si02 insulating layer. Compared to 
ohmic loss in the chip conductors, dielectric loss in the 
substrate can be neglected for the frequency and struc- 
ture size considered here poT |. as it is confirmed by our 
microwave simulation. 

We have checked that the potentials used in our simu- 
lation are robust against current and magnetic field fluc- 
tuations. We assume a relative stability on the level of 
10~ 5 for the static currents and fields. With available 
current sources and magnetic shielding, such a stability 
can be reached in experiments. The accuracy of the cur- 
rents (magnetic fields) specified in Sec. HVl is assumed to 
be better than 1 /iA (lmG). 

In order to study the robustness of the gate against 
noise on the time-dependent control parameters, we have 
simulated the gate with white noise on X(t) and uj±(t). 
The gate fidelity is reduced by ~ 10~ 4 for relative r.m.s. 
noise amplitudes n a < 10 -3 on the control parameters. 
For these values of n a , the small modulations of the ex- 
perimental parameters used for optimal control are well 
above the noise level, see Fig. |5] 



E. Two-photon transitions 

Another potential source of infidelities are two-photon 
transitions induced by the microwave. These can arise 
if more than one polarization component of the mi- 
crowave is present, as it is in general the case in mi- 
crowave near- fields. For tt 2 Ri A 2 , the two-photon 
coupling is characterized by an effective Rabi frequency 
^2 P h = R in iJ2 /(2A), where £l m , i <E {1,2}, are the 
single-photon Rabi frequencies of the transitions involved 
|37|. The detuning from two-photon resonance A2 P h = 
HbB{)/ (2K) is given by the Zeeman splitting due to the 
static magnetic field Bq in the trap center. For our trap 
parameters 



n 2 

max— i^- < 2 x 10~ 3 , 

A 2 — ' 
^2ph 



(33) 



and two-photon transitions are suppressed by the large 
two-photon detuning. 



G. Scattering length 

For an accurate theoretical description in one dimen- 
sion, the ID scattering length aip (fc) should be approx- 
imately independent of energy [38Ll39^ 



\aip(k) - aipjk = 0)| 
a 1D (k) 



< 1. 



(34) 



Otherwise, the phase shift will not be accurately pre- 
dicted by the ID calculation and a 3D calculation is nec- 
essary. In the experiment, the phase shift can never- 
theless be accurately measured and adjusted by slightly 
adjusting lu±. 

Due to the admixture of other hyperfine levels to the 
qubit states, the scattering lengths a°°, a° 4 , and a 11 are 
mixed. The order of magnitude of the changes in the 
qubit scattering lengths is 



Technical limitations 



'^-i&Ka^-aS 1 )-^- 0.02a 41 , (35) 



2A 2 



Since the conductor dimensions on our chip are smaller 
than the skin depth, the microwave current density in 



for 87 Rb, which typically gives Sa s ~ 10 4 a u , which is 
negligible. 
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VIII. CONCLUSION 

In conclusion, we have made a realistic proposal for a 
collisional phase gate using microwave potentials on an 
atom chip. The gate is implemented for a robust qubit 
state pair with experimentally demonstrated coherence 
and trap lifetimes r c ~ r t ~ 1 s at micron-distance from 
the chip. We have simulated and optimized the gate dy- 
namics for a chip layout which we specify in detail and 
which can be fabricated with today's technology. We 
found a gate fidelity of F = 0.996 at a gate operation 
time of T g = 1.1ms, taking many error sources into ac- 
count. With a total infidelity of the order of a few 10~ 3 , 
our gate meets the requirements for fault-tolerant quan- 
tum computation 0,0]. The gate fidelity of the present 
proposal is limited by wave packet distortion due to col- 
lisions and the dynamics of the atoms in a non-harmonic 
potential. We believe that these effects can be reduced 
by further optimization of the potential and better con- 
trol of the gate dynamics, at the expense of introducing 
more control parameters. The ultimate limit to the fi- 
delity will then be given by exp(— r g /r t ) ~ 0.999 for the 



chip layout discussed here. 

While the use of microwave potentials on atom chips is 
within reach, a major experimental challenge remaining 
is the deterministic preparation of single neutral atoms in 
the motional ground state of chip traps with very low oc- 
cupation probability of excited states. Proposals for sin- 
gle atom preparation have been put forward in [40L l4l| . 
An important prerequisite for single atom preparation is 
a single atom detector on the atom chip. The realization 
of such a detector is currently a subject of intense ex- 
perimental efforts |l6|, |42, |43j , with on-chip optical fiber 
cavities being a particularly promising system 0, 0] . 
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